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(N , 

Fisher-Hartwig asymptotics refers to the large n form of a class of Toeplitz determinants with 
' singular generating functions. This class of Toeplitz determinants occurs in the study of the 

spin-spin correlations for the two-dimensional Ising model, and the ground state density 
matrix of the impenetrable Bose gas, amongst other problems in mathematical physics. We 
give a new application of the original Fisher-Hartwig formula to the asymptotic decay of 
the Ising correlations above Tc, while the study of the Bose gas density matrix leads us to 
generalize the Fisher-Hartwig formula to the asymptotic form of random matrix averages 
over the classical groups and the Gaussian and Laguerre unitary matrix ensembles. Another 
, viewpoint of our generalizations is that they extend to Hankel determinants the Fishcr- 

, Hartwig asymptotic form known for Toeplitz determinants. 
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'Ph' Encomium 

> 

^ . birthday . 

^ \ Dyson's legendary works on random matrices are now standards in physics, mathematics and 

fields far(ther) afield. These works and, in particular his powerful log-Coulomb gas model, developed to 
liberate the mathematics where none yet exists, are the very essential tools in our ongoing pursuits 

1 Introduction 

Fisher-Hartwig asymptotics refers to the large n form of a class of Toeplitz determinants Dn{g\. By 
definition, the entries of the latter depend only on the difference of the row and column indices, and thus 

Dn[g\ = Aci[gj^k]],k=l,...,n (1-1) 

for some {gk\k=a,±i,±2,...- Crucial to the structure of the asymptotic form of (|l.lf) are analytic properties 
of the so called symbol 

oo 

g{e) (1.2) 



n— — OO 



^This radical appears in the works on Ramanujan 1 , who is very dear to Dyson. 
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or more particularly the decay of the Fourier coefficients of log^(^). Explicitly, let 

oo 

\ogg{e)= ^ Cpe^P'. (1.3) 

p—~oo 

Then if 

OO 

\p\cpc_p < oo (1.4) 

p— — OO 

a strong form of the Szeg5 Hmit theorem (see e.g. asserts that for n ^ oo 

oo 

Dn[g] = exp (nco + ^ kckC^k + o{l)j . (1.5) 
fc=i 

Two cases for which (|1.4ll wiU not hold are when g(9) has a jump discontinuity or a zero for some 
— TT < < TT. It is for such singular symbols (in the case of a zero it is the logarithm of the symbol which 
is singular) that Fisher and Hartwig [2] sought the asymptotic form of (|l.l|l . Symbols with singularities 
of this type have the functional form 

R R 
\ogg{e) = log a(6l)-i^6^arge'(^'-+''-^)+^a^log|2- 2 008(6* -61^)1 

r—l r—1 
R 

= loga(0) + (K + br) log(l + e*(«-(«-+-))) + (a, - K) log(l + e*(«'-+— »))) . (1.6) 
Here — tt < argz < tt and a{0) is assumed to be sufficiently smooth that if we write 

oo 

\oga{0)= Y Cpe'"' (1-7) 

P— — OC 

(cf. 11.3|l ) then the condition H1.4I) holds. By using data following from the fact that special cases of H1.6|) 
correspond to Toeplitz determinant expressions for the spin-spin correlation in the two-dimensional Ising 
model at criticality (see Section 2 below), the asymptotic form of which had previously been calculated 
0, Fisher and Hartwig ^6; conjectured that for some range of parameter values {ar}r=i....,R, {br}r=i....,R, 



Dn[g] e'=""e^ 
n — ^CJO 



where E is independent of n. Subsequently this was proved for various ranges of parameter values (see 
e.g. 0) and furthermore the constant was determined to be given by 



R 



g ^(a,+b,) log a_ (e,)g-(a,-6,) log a+ (S,) 



r=l 

(1 _ f^i{Ss-0r)y{a^+b^){as-bs) J| G(l + + br)G(l + Gr - K) 

where G is the Barnes G-function and 



^ •• - ' •• G(l + 2a^ 



oo —1 

loga+(0) :=^Cpe^''^ loga_(0) := ^ Cpe'^^ (1.10) 

p—l p—~oo 

Our interest is in applications and generalizations of the Fisher-Hartwig asymptotic formula (|1.8|l . We 
begin in Section 2 with an application of H1.8|l to the calculation of the asymptotic form of the spin-spin 
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correlation for the two-dimensional Ising model above criticality. In Section 3 the weU known equivalence 
of the Toeplitz determinant Hl.l|l to a random matrix average over the unitary group U{n) is revised. 
This average is in turn equivalent to the partition function of the one-component log-gas on a circle, 
subject to a one-body potential with Boltzmann factor g{9) at the special coupling (3 — 2. As such there 
is a natural generalization for couplings /3 > 0, and in the case &r = 0, r = 1, . . . , i? this can be used 
to predict the corresponding generalization of (|1.8|l . Moreover, in the special case a{9) = 1, i? = 1 the 
sought asymptotic form can be deduced from an exact formula valid for general ar,br. This can be used 
to extend the conjectured generalization of H1.8|) to non-zero br. 

In Section 4 we recall the problem of computing the asymptotic form of the density matrix for 
impenetrable bosons in Dirichlet and Neumann boundary conditions. This is immediately identifiable as 
an average over the classical groups Sp{N) and 0^{2N) respectively, with the function being averaged 
over having two zeros, and thus analogous to the random matrix formulation of the Toeplitz determinant 
(|l.l|l with symbol H1.6|l in the case R — 2, br = 0- We point out that the same class of averages over the 
groups 0^{2N + 1) or 0~(2N + 1) result from considering the density matrix for the impenetrable Bose 
gas in the case of mixed Dirichlet and Neumann boundary conditions. In |S] the sought asymptotics were 
calculated on the basis of a combination of analytic and log-gas arguments, and a Fisher-Hartwig type 
generalization (with b^ ~ 0) conjectured. The conjecture of can used to predict the asymptotic form in 
the case of mixed Dirichlet and Neumann boundary conditions. Moreover we show that this asymptotic 
form can be proved by making use of asymptotic formulas recently obtained 10 for Toeplitz -I- Hankel 
determinants 

det[aj_fe -f- aj+fe+i]j,fc=o,...,ra-i (l-H) 

in the case of singular generating functions H1.6|l . 

In addition to averages over the classical groups, the study of the density matrix for impenetrable 
bosons naturally leads to the question of obtaining the asymptotic form of averages over the eigenvalue 
probability density function for the GUE and LUE, in the case that the function being averaged over 
has zeros. Here the GUE denotes the Gaussian unitary ensemble of random Hermitian matrices, and the 
LUE denotes the Laguerre unitary ensemble of positive definite matrices with complex entries. These 
random matrix averages are equivalent to pure Hankel determinants 

/oo 
a{x)x" dfi{x) (1-12) 
-oo 

where dij.{x) — e^^^dx for the GUE and dii{x) — x"-e^^dx, x > for the LUE. Conjectures for such 
asymptotic forms are given in Section 5. The paper ends with some concluding remarks on the uni- 
versal form for Hankel asymptotics in Section 6, and attention is also drawn to the fluctuation formula 
perspective of our asymptotic results. 



2 Spin-spin correlations for the two-dimensional Ising model 

In the two-dimensional Ising model on a square lattice each site of the lattice exists in one of 

two possible states aij = ±1 with coupling between nearest neighbours in the horizontal and vertical 
directions. Explicitly, the joint probability density function for a particular configuration {aij} of the 
states on a {2N + 1) x {2N +1) lattice is given by 

^ N N~l N N-1 

P2N+l{{cr^J}) = C^p(^Ki (J,ja^+ij+K2 ^ ^ CT„cr^j + i) (2.1) 
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where Z2N+1 is the normaUzation. The spin-spin correlation function between the spin ctoo at the centre 
of the lattice, and the spin cr^.j. at site is, in the infinite lattice limit, defined as 



(cooo-i-j.) = lim <Jao(^i*j*P2N+i{{o-tj})- (2.2) 

N^OQ ^ — ^ 



Onsager knew of, but never published (see instead e.g. a Toeplitz determinant form for the case 
of H2.2|l for which (i*,j*) = {n,n) and thus lies on the diagonal. Explicitly 



(CTooCTnn) = det[a^ 



7— — l,...,n; 



where 

-1 + (l/fc)e-*''\i/2 



2 = -1 / h{B)e-'P^ de (2.3) 
27r J-^ 



/ i + i /fc e \ 

m i + (i/fc),.J ' fc-sinh2ifisinh2i^2. (2.4) 

Also, in the case of (|2.2|) with — (0,n) so that the two spins lie in the same row, Onsager and 

Kaufmann ^] expressed (|2.2|l as the sum of two Toeplitz determinants. A different approach to this 
problem was undertaken by Potts and Ward 13 , who obtained instead the single Toeplitz determinant 
form 

(ffooCTon) = det[a,_j],j=i,...,„, ^ j dO (2.5) 

where 

~ / (l + aie'^)(l + a2e-^) y/2 

''^^^ i(l + aie-«)(l + a2e^«)J ^^'^^ 

with 

ai := e~^^^ tanhifi, a2 



tanh Ki 

The formula obtain in 12 was shown to be identical to H2.3|l . (12.4(1 by MontroU, Potts and Ward |14j . 
We remark that if ai, a2 in (|2.t)|l and fc in (|2.4|) are regarded as parameters not specified by Ki, K2, then 
setting ai = 0, a2 = 1/fc in the former gives 12.4|l . 

A detailed study of the asymptotic form of H2.6|l was undertaken by Wu [7]. Indeed, it was the 
asymptotic form of 1)2. 5|l at the critical coupling 

ai < q;2 = 1 (2.7) 

obtained in j7j which, partially at least, inspired the formulation of the Fishcr-Hartwig asymptotic formula 
(|1.8|l "SI . To see how (|1.8|l relates to (|2.5f) with parameters H2.7|l , note 

logM^) =log(-i±^) +zarge-'^/2. (2.8) 



Q2=i \l + aie 

For |ai| < 1 this has the structure of 1)1.6(1 with 

/ 1 + ai p''^ \ 1/2 1 

«(^) = .J ' ^ = 1' &r = = 0, 6, = 

v 1 + aie / 2 

Recalling the definitions ()1.3)l (with g{6) replaced by a{6)) and ()1.1U)) . application of 1)1.8)1 implies 

l + aiNi/4 0FG2(l/2) 



*2 = 1 n^oo VI — Qfi/ 



where use has been made of the functional equation 

G{z + l)=T{z)G{z) 



(2.9) 
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in agreement with the resuh of Wu. 

The high temperature phase corresponds to couphngs 



ai < 1 < Q!2, oiici2 < 1- (2-10) 
In this case \ogh{9) is of the form with 

. ( ''+-!::;'; ^-"^-' 1 "^ « = , ... = -i. . o, «, ^ (2.11, 



(l + aie-»«)(l + e-'7a2) 

With R = l,br = — l,ar = we see that the Fisher-Hartwig asymptotic formula p.8|l breaks down 
because according to (|1.9|) the constant E contains the factor G(0) — and thus vanishes. To obtain 
the asymptotics in this case the approach taken in 7 was to relate it back to the original strong Szego 
theorem, multiplied by an auxilarly factor. Here we will show that by transforming (|2.fci|l . a form of 
\ogh{9) can be obtained which has the general structure H1.6|l but is distinct from the specification 
(|2.11|) . We will see that applying the Fisher-Hartwig formula then correctly reproduces the result of Wu 
for the leading asymptotic decay in the high temperature phase. 

For this purpose, let us introduce the notation f{d) = g{9) to mean that 



for some c independent of p. According to the definitions H2.6|l and (|2.11|) we have h{9) = e ^^a{9). Now, 
since with z = e**, h{d) is an analytic function of z in the annulus l/a2 < \z\ < a2, by Cauchy's theorem 



2tt J Jq 2mz 

for any simple closed contour encircling the origin in this annulus. Choosing C to be the circle with radius 
a2 (the outer boundary of the annulus) shows 



^' ~ U + (ai/a2)e-^V (I + e~'<' / alY'^ 

a2 (l + e-'7a2)i/2 Vl + (ai/a2)e-'«/ 1^1- 

This is of the form H1.6|l with 

11/ l + Q;iQ!9e*^ \ 1/2 3 1 

- —T. ' 2^^,A. , ' .J , i? = 1, 6, = a, = 9,. = -tt. (2.12) 

^2 (1 + e-'Va^)i/2 Vl + (ai/a2)e-*«/ ' 4' 4' ' ^ ^ 

Application of H1.8|l implies 



— n 

- aiy/\l - a,-^)-V4(i _ „^^^)-i/2 (2.13) 



in agreement with the result of Wu [7] . Moreover the Fisher-Hartwig formula 1)1. 8|l with R = 1 has been 
proved |15| for parameter values satisfying all three of the inequalities 

Reoi > 0, Re ai -t- Re 61 > — 1, Re ai — Re 61 > — 1. 

These inequalities are satisfied by the parameters in H2.12|l and so the Fisher-Hartwig formula provides a 
proof of (12331). 
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3 /^-generalization of the Fisher-Hartwig formula 

It is well known, and easy to verify, that the Toeplitz determinant 1)1. 2() can be written as a random 
matrix average according to 

n 

Dn[9] = {l[9{0i))^^^ . (3.1) 

\ / U(n) 

1 = 1 ^ ' 

Here U (n) refers to the eigenvalue probability density function for the unitary group 

n |e"^-e"^P, -n<e,<.. (3.2) 

^ ' ■ l<j<k<n 

As first noted by Dyson !T, H3.2|l is proportional to the Boltzmann factor for the one-component log- 
potential Coulomb gas on a circle, at the special coupling P — 2. From the log-gas viewpoint a natural 
generalization of 1)3.2(1 is the probability density function C/3E„ proportional to the Boltzmann factor for 
the same statistical mechanical system but with general coupling f3 > 0, 



(2^)"C„,^ ^^//^^J ' ' (r(/3/2 + !))»■ 

The identity H3.1|) then allows us to formulate a /3-generalization of the Toeplitz determinant p.2|) as the 
average 

n 

Di%]---{l[9i(^i))^^^ ■ (3.4) 

1=1 " 

Choosing g{9) according to in this we obtain a natural /3-generalization of the Toeplitz determinant 
with a Fisher-Hartwig symbol. In the case br = 0, r = 1, . . . ,R, the log-gas viewpoint can be used to 
conjecture the corresponding analogue of the asymptotic formula p.8|l . 
Let 

s'^^l'^. (3.5) 

l<j<fe<r 

As the first step, guided by both the log-gas viewpoint and the structure of ((1.8(1 . we conjecture the 
factorization 

From the work of Johansson plE], with the Fourier expansion of loga(0) specified by ((1.7(1 and 
assuming the coefficients satisfy (II .4(1 . it has been proved for general /3 > that 



(/3) 

-\\U\U }\ , , , ^, ^ ^. 



^co(n+Q) (2//3) Y,T=l f^f^kC-k 



Regarding the second ratio on the right hand side of ((3.6(1 . as first noted in ^|^] and revised in |2l)|. 
the log-gas viewpoint suggests that for n ^ oo we have the factorization 

TT |e»«. _ e'«. |^^.«^ " ^^^"^ ^ FT " ^^^^^^'^ (3.8) 

The large-n expansion of a ratio closely related to the product on the right hand side of 1(3.8(1 is known 
for (3 rational, in particular 

/3/2 = s/r, s and r relatively prime. (3.9) 
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Thus we have ^H] 
where 

ioji^(2g/r + i./r-p/5+l)G(^/r-p/s+l)- ^ 
Finally, the formula for C„.^ in H3.3|) together with Stirling's formula shows 

.(/3)r 



^f;-H^(r(/3/2 + l)r(n/3/2)-^'^/^ 



(3.12) 



Combining the above results gives the sought /3-generalization of the Fisher-Hartwig formula in the 
case br = 0- 

Conjecture 1. Let (3 be rational and of the form f^*. y\) . and let a{9) be as assumed for the validity of 
p. 71 ). For qjP > —1 we expect 



iV It 

Y\(ai9,)T\\e^'' -e^'^^l'i^P)) ^ e^o("+Ef.i 9.)„(/3/2) Ef.i (3.13) 

1=1 j=l ^ 

where, with Aq specified by 

R 

^ e"Sf=i«^^°sa(e3)g(2//3)i:r=i'=Cfc'=-fc Jl \e''^^ - e'^'\~^'^''^''W^Aq.. (3.14) 

l<j<k<R j=l 

It is of interest to extend Conjecture to include a factor 

in the average, and so obtain a /3-generalization of the Fisher-Hartwig formula for general parameters. 
Although we don't have a log-gas interpretation of the factor H3.15|l , the case R = 1 substituted in (|3.13|l 
with a{9) = gives an average which can be evaluated in closed form, and the corresponding asymptotics 
computed for rational (3. This together with the structure of the original Fisher-Hartwig formula 1)1. 8|l . 
p.9|l allows us to formulate the sought /3-generalization. 

Now, by rotational invariance, independent of the value of (j) 



(3) ■ (3.16) 

But, from the theory of the Selberg integral (see e.g. |3]), we know the right hand side of H3.16|l has the 
explicit gamma function evaluation 

fn{2cq,c) N TT fl/o ^Q17^ 

where /„(a,c) || , c := /3/2. (3.17) 



fn{c{q + b),c)fn{c{q-b),c) ' ' jj-^ {jc 
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For c € Z+ it was shown in Jl8j that 



Ua,c) ~ exp(anlogn)c""e-""n-(-i)"/^+"V^-n r^r^^ ''^w'^i^^ ' (3.18) 
"^oo J.X (^,((Q, _ + 1) 

while for r and s relatively prime 



G((a-p)/c+l)^ 



Using (pTTi^jl and IpTT^ in (pTTTjl it follows that for (3 rational of the form 

n 

TTe'/3fce,/2|^^g^e,|/3A _ (rn)^^^?'-''')^, ^ (3.20) 

/ CAE™ n^oo 

where 

Hio G(2g/r + ^/r-pAs + l)G(i./r-p/. + l) ' ^ 

Note that in the case 6 = this reduces to H3.11|l . (|3.10|) as it must. 

Knowing how, from the Fisher-Hartwig formula p.8|) . p.9|) . to generalize from the case i? = 1, general 
parameters, and the case general R but 6^ = (r = 1, . . . ,i?), to the case of general parameters and 
general R lets us use (|3.13f) and (|3.20|l to formulate a /3-generalization of the Fisher-Hartwig formula for 
general parameters. 

Conjecture 2. Let (3 be rational and of the form J^8. y\) . and let a{9) he as assumed for the validity of 
7| ). We expect, for some range of parameters {bj} and {qj}, 

n R 

where, with Aq^i, specified by 

R 

E(P) ^ g(2//5)i:~ ifcCfcC_fc -Qg-(g,+b,)loga_.(e,)g-(g,-b,)loga+(e,) 

r=l 

X II (i_e^(«-«^))-/3(9.+&.)(9.-f.)/2-Q^^^.,^^. (3.23) 

l<r/s<7? j=l 

4 Fisher-Hartwig asymptotics for averages over the orthogonal 
and symplectic groups 

A problem in mathematical physics which, along with the Ising correlations, motivated the Fisher-Hartwig 
formula p.8|) is the impenetrable Bose gas on a circle. If the circle has circumference length L, it was 
shown by Lenard |19| that the ground state density matrix p'^^^{x) has the Toeplitz determinant form 

Pn+i{x;0) = -^det[af_fe(x)]j-fc=i,...,Ar 

ap(x) := — / le^"^/-^ -f e*^||l-he*^|e-^'^d6l. (4.1) 
27r 
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The symbol in (|4.1|l is of the form (|1.6fl with 



a{0)^l, R^2, ai=a2 = -, 61=62 = 0, 6*1=0, 6*2=2^/1. 



(4.2) 



Now a fundamental issue relating to the Bose gas is the occupation Aq of the zero momentum state, which 
quantifies the phcnomcnum of Bosc-Einstein condensation (see e.g. j201). In the present system, which is 
translationally invariant, Aq is related to the density matrix by the simple formula 



Ao 



(4.3) 



For fixed < X < 1 one thus seeks the L 00 asymptotic form of p^^^{LX;0). In an unpublished 
work as of 1968, made available to the authors of 6 and subsequently published in 1972 Lenard 
obtained for the sought expansion 



p^+i(LX;0) -po 



y/2TT \Nsm{TTX) 



1/2 



(4.4) 



where po denotes the bulk density. Lenard obtained (|4.4|l as an upper bound, which soon after was shown 
to be attained by Widom [22] . Applying the Fisher-Hartwig formula H1.8(l with variables (j4.2|l reproduces 
g3I). 



According to H3.1|l the Toeplitz formula (|4.1() can equivalently be written as the U{N) average 



1 " 

^i(a;;0) = -(n 2sin 



1=1 



TTX til 

T " 2" 



2 sin ■ 



U(N) 



(4.5) 



The study of the ground state density matrices for the impenetrable Bose gas on a line of length L with 
Dirichlet or Neumann boundary conditions leads to formulas analogous to H4.5|l . only now the averages are 
with respect to the eigenvalue probability density functions for the classical groups Sp{N) and 0+(2A^) 
respectively (see e.g. |2 for the specification of these PDFs). Thus one has PU] 



PM+i{x]y) 



— sm — sm — Yy 2(cos cost^;) 2(cos cost/;) 



L \ 



1 ^ 



1=1 

TTX „ ^ 
COS ^ cost';) 2 (cos COS ( 



Sp(JV) 



0+(2Af) 



(4.6) 



Impenetrable bosons on the interval [0, L] with Dirichlet boundary conditions at a; = and Neumann 
boundary conditions at x — L also relate to a classical group. Thus from the fact that the single particle 
wave functions are given by 



2 TT(k-l/2)x 

Z"" L ' 



(fc = l,2,...) 



(the superscript M stands for "mixed" ) we see that the ground state wave function 

1 



C(xi, 



,xn) 



det[(^fc {x.j)]j^k=i,...,N 



has the product form 



,XN 



1 \N 



N 



Nl V2VI 



JJ^ 2 sin(7ra;i/2L) 2| costtx^/L — cosTrXj/Lj 

l<j<k<N 



1=1 
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The square of this quantity coincides with the eigenvalue PDF of the classical group 0^{2N + 1) with 
9 — TTx/L (for this we ignore the fixed eigenvalue at 6* = 0). From this fact, as in the derivation of H4.6|) 
detailed in .20., it follows that 

Pjv+ila;;?/) -^sm — sm— (^_[J_ |2(cos— - cos6',) 2 (cos cos 61,) (4-7) 



0+(2Af+l) 



Using a combination of analytic calculations based on the Selberg correlation integral |2H1 , and physi- 
cal arguments based on log-gas analogies, the large N , fixed x/L, y/L, N/L limit of the density matrices 
(|4.6|l was computed in |2j to be equal to 



D N „Gn3/2)(X(l-X))V8(y(i-y))i/8 



PN+l{x;y) - PN+l{x]y) P^= 1^ ^11/2 x = ri+co.,./w/2- (4-8) 



X = (l + coB xx/i)/2 
y = (l + coa,rH/i)/2 



Here we will show how recent rigorous asymptotic analysis [EI of Toeplitz + Hankel determinants 
with Fisher-Hartwig type symbols can be used to prove that pfj^^{x;y) exhibits the same asymptotic 
form (|4.8|l . We will also show how a result of 10 can be used to confirm the asymptotic form of a more 
general class of averages over 0+(2iV -I- 1) which can be deduced from a conjecture in |0, and how this 
conjecture in turn can be used to predict analogous asymptotics in the case of averages over Sp{N) and 
0+{2N). 

To begin we require a simple to verify identity noted in |23| . 
Lemma 1. Suppose g{9) = g{—9) and set gj — J"^ g{9)e'~^^^ d9. We have 

N N 

det[.9,-fe+.9,+fc+iUfe=o,...,^-i = {\{9{0i))^_ ={'Y{9{^-m)^^ ■ (4-9) 

Note that by the assumption on g{9) the matrix in (|4.9(l is symmetric. Also, the average in H4.7|l is 
an even function of 9i and corresponds to the special case 

g(9) = 2|cos— - cos^l 2| cos-^^ - cos6'| 

1j Lj 

= (|2 - 2 cos(6l - — )| |2 - 2 cos(6l + — )| |2 - 2 cos(6' - ^)| |2 - 2 cos(6l + — )|) (4.10) 
V L L L L J 

of (|4.9|l . We observe that (|4.10|) is an example of a symbol of the form (|l.t)|) . Fortunately, recent rigorous 
works [301 [TU) have determined the asymptotic form of the Hankel -I- Toeplitz determinant in H4.9|l for all 
symbols (|1.6|l . with the restriction that for g{9) even (the case of interest in relation to (|4.9|) ). 9r ^ 0, ±7r. 
Let us recall the result of fOl Thm. 6.1], simplified so that it relates to the even case of (ll.6|l with each 
hr ^ 0. 

Theorem 1. Let 

R 

logg{9) =loga(6i) + ^a^(^log(2|cos6'-cos6'^|)), (4.11) 

where a{9) is an even periodic function with the property that the Fourier expansion of its logarithm jl.T)) 
satisfies j j.^l ), together with some technical assumptions (for the latter, which may not be necessary, see 
IWH). We have 

det[g,,k + gj+k+ihk=o,...,N^i - e(^+Sf-i'^^)^°(27V)S?-i'^'i? (4.12) 
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where 



^\ G{l + 2ar) 



xIl M+jLHit..,., n (|l-e^^^^-=)||l-e^(^-«-=)|)-^"'-"\ (4.13) 

\<T<S<R 



Recalling (|4.9(l it follows from Theorem that 

AT 

n|2(cos — — cos6',-)l |2(cos — — cos 6',) 
' ^ L ^ L ^'-l o-(2N^\) 

M , ma:/L|l/4 M I 7rii//L|l/4 -i 



Substituting this in H4.7|l shows 



(27V)l/2 4 |2 _ g27ri:x/L|l/4|j^ _ g27ria/L|l/4 



g4(3/2) (x(i - x))i/8(r(i _ r))i/8 



/2]v |x-y|i/2 



x.a,oo= (4-15) 

Y' = (l + coa,ry/t)/2 



thus rigorously establishing the asymptotic form (|4.8|l derived, but not rigorously proved, in ^ for the 
cases of Dirichlet and Neumann boundary conditions. 

The eigenvalue distributions for Sp{N), 0+(27V), {2N + 1) and 0+{2N + 1) are proportional to 

AT 

J|(l + cos6i,)^i(l-cos6'0^' n {cos9k-cosejf, < Oi < n (4.16) 

1=1 l<j<k<N 

for (Ai, A2) — (1, 1), (0, 0), (1, 0) and (0, 1) respectively (our convention is not to include the delta function 
corresponding to a fixed eigenvalue, nor the delta functions corresponding to the conjugate eigenvalues). 
Thus to obtain the asymptotics of the averages in (|4.t)|l it is sufficient to obtain the asymptotics of 

N R 2 

(nn(2|cos^.-cos0.|) ^) (4.17) 

1=1 r=l 

where CAr(Ai, A2) refers to the normalized form of (|4.16|l . In the cases (Ai, A2) ~ (1, 0) or (0, 1), due to the 
identity (|4.9|l . we can read off the asymptotic form from H4.12|l . But for general (Ai,A2) the asymptotic 
form of (|4.17() is not included in Theorem ^ Instead, we will use a conjecture from P| to formulate the 
result. 

Let us first recall the conjectured asymptotic form from |9|. Define 

HnMM[f{^)]--= dxi--- dxnY[f{xi)xf'{l~xi)^^ l[ \xk-Xj\\ (4.18) 

•^0 -^0 1=1 l<j<k<n 
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Then the argument given in |9j predicts^ 



KEtl9r + (Al+A2)/2 h{x) 



exp 



Wl-a;)]V2 



dx 



R 



I 2^['''lr + q;) log 2n 

r=l 



K 



(4.19) 



where 



K = 



j=l l<j<k<R 



X exp 



[47r2 



(x(l - a:))V2 



1 G'(gr + 1) 
,i-\^'-G(29. + l)- 

In preparation for relating this to the average (|4.17|) let 



h 



i nO. 



We then have that 



(^i(l + cos0)^ = Co + 2 Cn cos ) 

n— 1 

1 1 °° 

-(Ai/l(0) + \2h{l)) = -(Ai + A2)co + Cn(Ai + (-1)"A2) 

n—1 

while H4.21|l together with the cosine expansion 

oo 2 

log(2 1 cos ~ cos (t)\) — ~ — cos nO cos nt/) 



shows 



Also, as noted in [3], 



47r2 7o (x(l - x))i/2 



Mx) /^'(y)(y(i-y))^/^ _ ig^, 



X - 1/ 



G(n + 1 + a)G(n + 1 + b)G(n + 1 + a + b) 



G(l + a)G(l + 6)G(2n + 1 + a + 6) 



Since [ID 



log 



G{n + 1 + a) 
G(n + 1 + 5) n- 



(& - a)n + ^— log(27r) + ( (a - 6)n + 



a2-62 



logn + o(l) 



(4.20) 

(4.21) 

(4.22) 
(4.23) 



(4.24) 



(4.25) 



(4.26) 



^Unfortunately there are a number of inaccuracies in the reporting of the conjecture in 0. The term 
H^^Y^R denominator on the left hand side of (14.1911 has mistakenly been written as -ff„,Ai,A2[l] 

in equations (90), (94), (96) and (97); the factors IljLi 2/^ "^^'^ (1 — yj)~^'^'^' are missing and should be paired 



with W 



l<j<k<R 



~ yj\ ^'^"^ throughout; and the term e 



-(Ai+A2)[h(0) + h(l)l/4 



in (96) and (97) should read 



-(Aih(0) + A2h(l))/2 
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we deduce 

IT ril ')4nQ+2Q^+2Q{a+b) 

^^•^•bW ^ t (A 27) 

i/„+Q,,,6[l] {27Tn)Q • ^ ■ 

Finally we note that under the change of variables 

XI = i(l + cos6l() 

the integrand in l|4.18|l contains as a factor the (unnormalized) eigenvalue probability density function 
(jllBI) . Explicitly, with h{9) h{^{l + cos 6*)) we have 

H.^..M[e^y^\yr-.\^^^] =(neM«')n|i(cos0.-cos.O^n . (4.28) 



Cn (Ai ,A2) 



Making use of 14.22|I - H4.28|I shows that the asymptotic formula ()4.19|l for the integral (|4.18|) is equiv- 
alent to an asymptotic formula generalizing Theorem ^ 

Conjecture 3. Let \oga{9) have the Fourier expansion jJ.Tl ), with coefficients satisfying jl-j}) - W^e 
expect that for N —i- oo 

N R 

(T\a{ei)T\\2{cos(j)r-cosei)\^''A _ e(^+^-i"'-)'=°(2^)^^=^'''^ (4-29) 

\ ^ ^ ' Cjv(Ai + l/2,A2 + l/2) 

where, with E specified by j^--?^ , 
R 

We can apply some checks to (|4.19|) . As already remarked, with (Ai,A2) = (1,0) the probability 
density function Cn{M, A2) coincides with the eigenvalue probability density function for 0~{2N + 1), 
and (|4.29|l must coincide with H4.12|l . as indeed it does. Also, changing variables 9i n ~ 9i and 
interchanging Ai and A2 leaves (|4.16|l invariant, and thus the average H4.17(l invariant if we also put 
(/)^ h- > TT — (fjr, a{9) a{TT — 9) (and thus c„ ^ (— l)"c„). Recalling the definition (|4.13|) of E we see that 
1)4. 30(1 exhibits this symmetry. Another check follows from a factorization identity, relating an average 
over the unitary group to a product of averages over the orthogonal and symplectic groups |25l 1231 1^ . 

Proposition 1. With g{9) = g{—9) we have 

2N+1 N+l N 

{ n ."C')),,,.,,,, - { n »w)>o.,»„,(n^'(''.)),„„, c-'d 

With 

R 

3(6*) = 0(6*) J|(2| cos6' - cos(/)^|)2°- 

in (|4.31(l we see that the conjectured asymptotic form 14.29|l for the right hand side is consistent with 
the Fishcr-Hartwig formula (|1.8|) . 

The identity H4.31() is also of interest for providing an exact formula for the product of the density 
matrix in Dirichlet boundary conditions and in Neumann boundary conditions. Thus recalling H4.6|l we 
see from (|4.31ll that 



2N+1 

— sm — sm — ( II 2 (cos— -cose,) 2(cos — - cos 6I/) ) = pj^^^{x,y)pj^^^(x,y). 



1 . TTX . nil I "-f-T L , TTX „ ^ Tra; 



1=1 
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5 Impenetrable bosons in a harmonic trap and random matrix 
averages over the GUE and LUE 

From a physical viewpoint the most relevant setting for the impenetrable Bose gas is confinement by 
a harmonic potential (see |2(J[ I27| and references therein) . Then the ground state wave function V'o^ is 
proportional to 

N 

1=1 l<j<k<N 

and IV'o^P is identical to the eigenvalue probability density function for the Gaussian unitary ensemble 
of complex Hermitian matrices. The combination of log-gas arguments and analytic calculation based on 
the Selberg correlation integral used to analyze H4.18|l was used in [27] to analyze the asymptotic form of 

e ^^=1 7- (5.1) 

where 

/OO /-OO ^ 

dxi--- dx^ []/(x,)e-"'^' n \^k-xj\^ (5.2) 

in the special case R = 2, qi = q2 — 1/2 which specifies the ground state density matrix. As our first 
point of interest we will generalize this calculation to general R and Qr {qr > —1/2). 
The log-gas perspective suggests the factorization 

(5.3) 

Next, from the theory of Selberg correlation integrals |2H1, for G ^>o we have the duality formula [221 

CTV.ail] G2q^,a[l] 

It is a fairly straight forward exercise, detailed in |29[ section 5.3] for a related problem, to use the saddle 
point method to compute the large N expansion of the integral on the right hand side of 1)5. 4|) . One finds 

e-''"'^^^^G,^^V2N{iyr + ^^r]- (^''){G,^A^]fe-^^^2-'^^^m)-^Hl-y'y^'- (5-5) 

N-^Qo \ q^ / 

Also, we know (see e.g. |3j) the exact evaluation 

G„,i[l] = 2-"("-i)/27r"/2G(„ + 2), (5.6) 
which together with the asymptotic expansion 14.26|) implies 

\lN]Sl-q^ (5.7) 



Combining H5.4|l - 15.7|l gives the asymptotic formula 

2Na.yl GN,^[\^~Vr?'^^] ^ (g. + 1) 2^^^^^ g,, 

G^+,.,vw[l] ^--G(2g. + 1) TT.^ ^'-'^ 

which we have proved for qr G Z>o, and conjecture as being valid for all q^ > —1/2. This same result can 
also be deduced from results in 31 . Substituting this in 1)5. 3|l gives the sought asymptotic form of (|5.1)l . 
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Conjecture 4. Let qr > —1/2, and let GN,a[f] be given by \5.SJ^ . We expect 



G 



n 



\yj -yk\ 



G2(g, + 1) 229^ 



l<j<k<R 



nG(2g, + l) TT?^ 



(5.9) 



An extension of H5.9|l can also be formulated. Let a{x) be analytic on [—1, 1]. Then it has rigorously 
been proved that |TB] 



G 



N,V2N 



[ga(rz:)j 



G 



N,V2N 



[1] 



-1 a;-?; / 



The structure of H4.19|l in the case Ai = A2 = suggests how (|5.1Uf) can be combined with l|5.9|l to 
generalize the latter. 

Conjecture 5. Let a{x) be analytic on [—1, 1]. It is expected that 



G 



[1] 



~ RHS ilTTTIt 



(5.11) 



We remark that in the special case a{x) = kx, Conjecture |31 can be reduced to Conjecture 01 To see 
this, use completion of squares to note 



G 



N.V2N 



kx 



n 



x-yr 



ri^G 



N.V2N 



[n 



X H Vr 

2N ^ 



According to Conjecture 01 we have 



-Y.^^i2Nq^iy^-k/4Ny 



-y P-Jr 



and thus 



G 



e*^'/*^ RHS lirrm 



Since H4.24|l gives 



1 



dx ■ 



aix) r ^^a'(y)(l-y2)i/2 



a(a:)— /ca: 8 



47r2 7_i (1 - a;2)i/2 - a; - y 

this is in agreement with H5.11|l . 

Let us now turn our attention to a variation of the impenetrable Bose gas in a harmonic well, which 
also has the features of being related to a random matrix ensemble. In reduced units, the Hamiltonian 
for the system is 



H 



i:|.+i:(«v-i)^ 



N 



, x'^ 

2 ' -^J 



Xj > 0. 



(5.12) 



Thus in addition to the harmonic well, the particles are restricted to the half line by a repulsive potential 
(requiring a' > 1) at the origin proportional to l/r^. This is the non- interacting case of the so called 
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type B Calogero-Sutherland Hamiltonian 32 , for which the interacting case has pair repulsion. The 
ground state wave function for H5.12|l is proportional to 



j^^-.?/2(^2).72 Y[ \xl~x]\. (5.13) 

l=\ l<j<k<N 

We recognize the square of the ground state wave function as being identical to the probability density 
function for the singular values of n x N complex Gaussian matrices with a' = n — N + 1/2 (see e.g. 
Ch. 2]. Changing variables xf i-^ xi this is refered to as the Laguerre unitary ensemble. The problem of 
computing the asymptotic form of the density matrix for this system suggests analyzing the asymptotic 
form of the more general quantity 

where 

LN[f]--= dx,--- dx^, n/(^')^^e-4'^'' n l^fe-^'l'' (5-15) 

■^0 ''^ 1=1 l<j<k<N 

in analogy with H5.1|) . 

From a log-gas perspective, the integrand in (|5.15|) corresponds to a one-component system interacting 
on the half line x > 0, subject to a one-body confining potential 2iVx^ — (a' — 1/2) log x. In addition to the 
electrostatic energy — log jx — x'| at the point x due to the interaction with a charge at a;', there is also a 
term — log |a; 4- x'| due to the interaction with an image charge at ~x' (outside the system, since x' > 0). 
In keeping with the image charge interpretation, for each charge at x one requires a term — ^ log \ 2x\ to 
account for the interaction between a charge and its own image (the factor of 1/2 is because this energy 
is shared between the charge and its image, the latter being outside the system). From this viewpoint 
we can interpret 

l<3<k<R r=l ^^''='^N+Z?^,qA^\ 

as a ratio of partition functions for log-gas systems, and analogous to H5.3|l we expect the factorization 
into 



n(2y,)^?e--4^^=?|,.|(-'-i)-M^-^ (5.16) 



for N oo. 

To analyze (|5.16|l in the limit N oo we make the change of variables xf i— > xi and introduce 

/■oo /"OO ^ ^ 

LN.clf]-^ dxi--- dxjv n/(2^0a^r"'^'e-^"' H i^k - x,)^ (5.17) 



1=1 l<j<k<N 



SO that it reads 



n(2;/.)^S^e-^-^^^^y^''^)^- ^^'.^^['^ ~ y'^'"^^ . (5.18) 
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To proceed further, we use the fact that for q £ Z>o we have the duahty formula |33 | 



i7V,c[l] 



M2,(a,7V) 

at-^a-\-2q 
1/2 nl/2 2g 



-1/2 J-1/2 77 ,^~-T:^ 



|g27ria:fc _ ^2TTiXj |2 



"1/2 "'-1/2 l<j<k<2q 

(5.19) 

where on the right hand side a = a' — 1/2 and 

r.l/2 /.1/2 



\^2TTiXk _ ^2ivixj 1 2 



-1/2 J -1/2 ,</<t 

-G{n + 2) (5.20) 



"1/2 "'-1/2 l<j<*:<n 

G{n + l + a + b) G{l + a) G{l + b) 



G{l + a + b) G{n + l + a)G{n + l + b) 

(for the last equality see e.g. [Sj). If we suppose temporarily that a e Z>q, the right hand side of (|5.19|) 
with c = AN can be written as the contour integral 

where C is any simple closed contour which encircles the origin. To analyze this for iV — > cx), following \,VA\ 
where the case 9=1 was considered, we note the TV-dependent terms in the integrand have a stationary 
point when 

z^z^:=-\±r\{\/t-lY/\ (5.22) 

By deforming the contour C to pass through z+ for q of the integrations, and to pass through z_ for the 
remaining q integrations, we readily deduce from the representation (|5.21|) of H5.19(l that 

e-,.4^.,^y(2a'-l),.(22;,)^i!L^,4^[|x - y^^..] 



^/2q,\ ^-4NyUA^^+z^)+Nq. log |l + l/.+ |^ - z_ 1 1 + ,^\1r(2a'-l) f J_\ 

\ qr ^ k+l V27r/ 



/ 1 1 



2\zl [l + z+Y 
Now, with t^yl'in 



''\G,Al]?. (5.23) 



z+ 

+ (5.24, 

SO the right hand side of 1)5.2311 simplifies to 

Furthermore we know (see e.g. |H]) 

~ _ ^_^^-^(.-_i/2) g(^ + 2)G(a' + iV+l/2) 
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and making use too of (|5.22(l it follows from the asymptotic expansion (|4.26|l that 

iw+9.,4iv[l]M2,,(a'-l/2,iV) 

Substituting H5.27|l in H5.25|l . evaluating Gq,^!] therein according to H5.6|l and simplifying we obtain the 
N oo expansion 

proved for qr € Z>o and expected to be true for all qr > —1/2. Substituting this in (|5.1t)|l gives, as a 
conjecture, the sought asymptotic form of H5.14|l . 

Conjecture 6. For N oo, and assuming qr > —1/2 for each r = 1, . . . , i?, 



n. 



■ ^ " r— 1 V / 



l<j<k<R 



It is of interest to extend (|5.29|) in an analogous way to how (|5.11|l extends ()5.9|l . First we use 1)4.19(1 
with Qr = 0, and (|5.1U|) to conjecture that for a(x) analytic on [0, 1] 

^^^exp(^(^+(^"'-^)/^) raix)VT^^dx) 

^».(^r-^/>.^^^^t^)--'' 

Combining this with (|5.29|l as in (|5.11(l gives us the LUE analogue of Conjecture 
Conjecture 7. Let a(x) be analytic on [0, 1]. It is expected that 



r=l 



- (rRS^EM) fRHSdHa^e"^ (5.31) 

We can check the consistency of ((5.11(1 and ((5.31(1 . For this we make use of a factorization identity 
analogous to Proposition ITI 

Proposition 2. Let g{6) = g{—9). We have 



where 



/OO pOC ^ 

dx,--- dx^l[g{xi)\xi\Pe~^'^' JJ (4 ^ ^ 



(5.32) 



1=1 l<j<k<N 
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Let a{x) be even and choose 

R 

r=l 

According to Conjecture 



TV— ►oo 



,4 ^ /•! 

exp 



X n Is? -»ir*"» rid -»;)'•" i2»r"'= 

i<j<k<R r=l 
R 

X 



r— 1 ' 

For the right hand side of (|5.32|) as imphed by Conjecture [7| to be consistent with this we require 

o4y:5_, G2N,V4nW '^2(Af+i:«^,g,0,v/4]v[l] 



But the method of derivation of (|5.32|l given in |21] shows that for general 



Lr.\a'=0Ln\a' = l (n!)^ (™)l/2' 

verifying (|5.34|) . 

Let us now apply Conjecture[7|to the calculation of the density matrix pj^^^{x, y) for the state (|5.13|) 
with + 1 particles, 

' {x,y):^^e-'^'/^-y''\xvf 



AT 

„2 2^2 

l=\ l<j<k<N 



poo poo 

/ dx,--- dxr,Y[xfe-^'\x'-xf\\y'-yf\ J] ( 



xl - x-]f (5.35) 



where 

l-oo roo N+1 



N+1 ■= dxi-- ■ dxN+i Y\_ '3"^' Yl_ i^k - 

•^0 "^0 1=1 l<j<k<N+l 

In terms of the average H5.15|l we thus have 



On the right hand side we can apply Conjecture[3with R = 2, qi = q2 = 1/2 and so obtain the asymptotic 
form 

2^/iVp^+l(2ViVX,2^/iVr) ^ 2V7V^^^^^^^^(1 _ x^^^il - Y^Y'\ (5.36) 
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The asymptotic form (|5.36() can in turn be used to specify the occupations Aj of the low-lying effective 
single particle states which by definition satisfy the eigenvalue equation 

PN{x,y)(j>j{y)dy = Xj(j)j{x). (5.37) 

Thus, with X = 2\/NX, y — and j fixed, introducing the scaled effective single particle states 

El EDI 

substituting H5.38|l and using the fact that pj(j{x, y) is supported on a:, y € [0, we obtain the explicit 

integral equation 



2 

where 



A^. = ViV^^^A,. (5.39) 

TT 

We see immediately that the occupations of the low-lying effective single particle states are proportional 
to ^/N , as has been found for the impenetrable Bose gas in periodic boundary conditions ^lEOl, in a 
harmonic trap |34[ I27| and in Dirichlet and Neumann boundary conditions (H]. An appropriate analysis 
similar to that undertaken in [271 Appendix B] gives the same upper bound on Aq as found for the same 
quantity in the case of the harmonic trap [23, but a detailed analysis of (|5.38|) remains. 



6 Concluding remarks 

6.1 Universal form for Hankel asymptotics 

Analogous to H3.1|l . Hankel determinants are related to log-gas partition functions according to the formula 

det [aj+k\j,k=Q,. ..,n~i 

-1 />oo />oo ^ 

1=1 l<j<k<n 

=: A„(e-"^(-))[a(x)] (6.1) 

where 

ap= a{x)xPe-"^'^''Ux. 



For V{x) an even degree polynomial independent of n with positive leading coefficient and no real zeros, 
it was proved by Johansson pi 6) that 

yl„(e-"^(^))[e°(^)] 
A„(e-"^(-))[l] 

exp (n I a{x)p{x) dx\ exp [ dx ^^^] [ dy*^-^ — ^— — — —\ . 

"-'"^ ^ Jci ^ ^47r^7ci J [x ~ ci){c2 - x) Jc, x-y ) 



(6.2) 



Here p(x) is the scaled density in the log-gas system corresponding to An{e "^'^^))[1], supported on 
[ci , C2] and normalized so that 

p(x) dx = \. 
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The asymptotic formula H5.10|l corresponds to the special case V{x) = \x^, p{x) = fVl-2:2 of dOi . 
To extend Conjectm'e[Slto more general V this suggests we simply write the latter in terms of p{x). 

Conjecture 8. Under the conditions of the validity of 1^6. we expect 



l<]<k<R r=l ^ ' 

We remark that in the case R = 1, e"^^^ = 1, this conjecture (together with some corroborative 
analysis) was formulated earlier by Brezin and Hikami |H3 (see also ^^). 

Conjecture dean similarly be extended, although we work with the quantity (|5.17(l in favour of (|5.15(l 
so as to have a Hankel determinant interpretation according to ()6.1|l . In the log-gas system corresponding 
to (|5.17() one has p{x) = ^^^j^ (1 — x)^!'^ . Recalling the equality between (|5.16|l and H5.18|l . and writing 
?/r ^ Vrj 0'[x^^'^) i~> a{x) we see that Conjecture can be rewritten to imply 



n 



ia' -1/2),^ ~iN,^y^ LjV^Af [e°(^) Ilf^l ^ " yA^""] 



e n - y^-''^'' n rjo^' lu C^^Ny^-^'ipiyrW'-, (6.4) 

l<j<k<R 



N^oc _ J-J- G(2gr + 1) 

r— 1 ^ ^ 



thus assuming the universal form (|6.3|l and suggesting the following analogue of (|6.2|l and Conjecture |H| 

Conjecture 9. Let V{x) be a polynomial independent of n, with positive leading coefficient and no real 
zeros on [0, oo). Let 

-| /"OO /"OO ^ 

i„(x"e-"^(^))[a(x)] - / dxi x?e-"^(^^) • • • / dx„ <e-"^("") J] «(^') 11 (^fe-^j)'- 

(6.5) 

Analogous to 1^6. ^) we expect that 

i„(x"e-"^(^))[e''(^)] 
i„(x"e-"^(^))[l] 

~ exp (n [ a{x)p{x) dx \ exp ( — - [ dx — ^=LJ^= [ dy — (^)v^^( ^ (6.6) 
«-.oo V Jo / V47r^ 7o ^a;(c2-x)Jci - 2/ ^ 

where p{x) is the scaled density in the log-gas corresponding to A„(a;"e^"^^^^)[l], with support on [0, C2]. 
Furthermore, with the same meaning of p{x), we expect 

As a final comment on this point, we note that the universal form given by the right hand side of 
1)6. 3|l is also exhibited by the Fisher-Hartwig formula (|1.8|) . Thus, with Zr :— e'^'' we see that 



p(y)=N/2„ 
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6.2 Further Toeplitz + Hankel structures 

The identity (|4.9|) of Lemmas has counterparts for averages over Sp{N) and 0+(2A^) |23j. 

Lemma 2. Suppose g{6) = g{—6), set gj = J^^ g{9)e~'^^^ d9, and let C^i^i, ^2) refer to the normal- 
ized form of i4.1S\) . We have 

N N 

det[a,_.+a,+.U=o....,^-i = (Uai^^^^) o.^^n) = i^^'^^^^) c.io,o) 

N N 

det[aj-k- aj+k+2]j,k=o,...,N-i = ( H ^^^^'Ve ( 11 ■^(^j)/^ ,^■ (6-^) 

^ , ' Sp{N) ^ , / Civ(14) 

Choosing g(0) as in (|4.11|l . Conjecture O gives the asymptotic behaviour of the right hand sides in 
Ht).8|l . and thus the conjectured form of these Toeplitz + Hankel structures. 

6.3 Fluctuation formula perspective and future directions 



Let p := p{xi, . . . jX^r) be an A^-dimensional probability density function. The stochastic quantity A 
defined by 



SiLi with the {xj} sampled from p, is referred to as a linear statistic. Its distribution -Pa(^) is 



N 



PA{t)^(^6{t-J2aix,))) , (6.9) 



and taking the Fourier transform of this gives 

N 



PA(k) = m e*'='^(^^M . (6.10) 



0=1 



The structure of the average (|6.1U|) is common to the averages studied in this paper. As an illustration 
of the content of the asymptotic formulas from this viewpoint, consider Johansson's result Ij^.Tfl . Written 
in terms of the average (jXT|) with g{e) = e'^<^\ it reads 

£)(/3)jg»fca(e)] ^ gifecong-(2//3)fc2 2~=i "c„c_„ (6.11) 
n — >oo 

where {c„}„=o,±i,... are the Fourier coefficients in the expansion of a[6), 

oo 



n— — oo 



A key feature of the exponents in the exponentials on the right hand side of H6.11|l is that they form a 
quadratic polynomial in k. Thus substituting this in (|6.9|l and taking the inverse transform gives the 
Gaussian distribution 

''M,^::^^^.-"-'"-"'- (6.13) 



with 



4 °° 



fi = nco, a'^ = -^pcpc-p. (6.14) 

P p=i 

As noted by Johansson (5j, in the case f3 — 2 this gives a Gaussian fluctuation formula interpetation 
of Szego's theorem. A perculiar feature is that although the mean is proportional to n, the variance is 
0(1), so fluctuations are strongly suppressed. It is formulas of the type H6.13|l . (|6.14|l which led to the 
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successful theoretical explanation of the phenomenon! of universal conductance fluctuations in mesoscopic 
wires (see e.g. \dG\). in which the conductance — an order N quantity — is written as a linear statistic 
of certain eigenvalues and is shown to have 0(1) fluctuations with variance given by an analytic formula 
of the type H6.14|l . 

All our generalizations of the Fisher-Hartwig formula involve a term of the form e*^ as the first 
correction to the leading order behaviour e"^"". However again when written as an average of the type 
Ht).10|) the exponential of a quadratic in k again results. Consider for example H3.13|l . With {c„} specified 
by (|6.12|) we have 

/^ikai9)+ikl3Y.j = i log le^'-e'^j | \ ^ gifcco«g-fe2(/3/2)(5:f^, qf) log n 

\ / C/3E„ II-+00 

and thus, as first noted in [37], with 

N R 

A = ^(a(0O + /3E9^1og|e^''-e^^^l) 

the asymptotic form of the corresponding distribution is given by the Gaussian ()6.13fl with 

R 

fi = ncQ, = /5( X! ^l) n. 

Thus the variance diverges logarithmically. This class of Gaussian fluctuation theorem has found use 
in the application of random matrix theory to the study of the statistical properties of the zeros of the 
Riemann zeta function |38l 15^ . The study of the statistical properties of the zeros of families of L- 
functions requires averages over the different classical groups |4(J[ 1411 1421 021 ■ We might anticipate that 
our new results of Section 4 will find application in this topic. 

Of course it remains to prove the conjectures of this paper. Of these, ConjectureHlis the most general, 
as it involves Fisher-Hartwig type parameters {g^}, {bj} as well as the log-gas type parameter /3. It is 
also of interest to extend Conjectures|31|Hland|51to this level of generality. Another direction of generality 
is to extend the domain of integration from a circle or line to a two-dimensional region [441 145| . 

Acknowledgement 

The financial support of the Australian Research Council is acknowledged. 

References 

[1] B.C. Berndt, H.H. Chan and L-C. Zhang. Amer. Math. Monthly 104, 905 (1997); Acta Arithmetica 
87, 145 (1998) 

[2] F.J. Dyson. J. Math. Phys. 3, 140, 157, 166, 1191, 1199 (1962); 13, 90 (1972); Commun. Math. 
Phys. 19, 235 (1970); 47, 171 (1976) 

[3] P.J. Forrester. Log-gases and random matrices, www.ms.unimelb.edu.au/~matpjf/matpjf.html 

[4] H. Widom. Adv. Math. 21, 1 (1976) 

[5] K. Johansson. Bull. Sci. Math. 112, 257 (1988) 

[6] M.E. Fisher and R.E. Hartwig. Adv. Chem. Phys. 15, 333 (1968) 



23 



[7] T.T. Wu. Phys. Rev. 149, 380 (1966) 

[8] E.L. Basor and K.E. Morrison. Linear Algebra Appl. 202, 129 (1994) 
[9] P.J. Forrester, N.E. Frankel and T.M. Garoni. J. Math. Phys. 44, 4157 (2003) 
[10] E.L. Basor and T. Ehrhardt. Oper. Th. Adv. and Appl. 135, 61 (2002) 

[11] B. McCoy and T.T. Wu. The two-dimensional Ising model. Harvard University Press, Harvard, 1973 

[12] B. Kaufman and L. Onsager. Phys. Rev. 76, 1244 (1949) 

[13] R.B. Potts and J.C. Ward. Prog. Theor. Phys. (Kyoto) 13, 38 (1955) 

[14] E.W. MontroU, R.B. Potts and J.C. Ward. J. Math. Phys. 4, 308 (1963) 

[15] A. Bottcher and B. Silbermann. Math. Nachr. 127, 95 (1986) 

[16] K. Johansson Duke Math. J. 91, 151 (1998) 

[17] P.J. Forrester and C. Pisani. Nucl. Phys. B 374, 720 (1992) 

[18] P.J. Forrester. Phys. Lett. A 163, 121 (1992) 

[19] A. Lenard. J. Math. Phys. 5, 930 (1964) 

[20] P.J. Forrester, N.E. Frankel, T.M. Garoni and N.S. Witte. Commun. Math. Phys. 238, 257 (2003) 

[21] A. Lenard. Pacific J. Math. 42, 137 (1972) 

[22] H. Widom. Amer. J. Math. 95, 333 (1973) 

[23] J. Baik and E.M. Rains. Duke Math. J. 109, 205 (2001) 

[24] E.W. Barnes. Quart. J. Pure Appl. Math. (Germany) 31, 264 (1900) 

[25] H.S. Wilf. J. Comh. Theory Ser. A 60, 155 (1992) 

[26] P.J. Forrester. Evenness symmetry and inter- relationships between gap probabilities in random 
matrix theory. Preprint (2003) 

[27] P.J. Forrester, N.E. Frankel, T.M. Garoni and N.S. Witte. Phys. Rev. A 67, 043607 (2003) 

[28] J. Kaneko. SIAM J. Math. Anal. 24, 1086 (1993) 

[29] T.H. Baker and P.J. Forrester. Commun. Math. Phys. 188, 1371 (1997) 

[30] E.L. Basor and T. Ehrhardt. Math. Nachr. 228, 5 (2001) 

[31] E. Brczin and S. Hikami. Commun. Math. Phys. 214, 111 (2000) 

[32] T. Yamamoto. J. Phys. Soc. Japan 63, 1223 (1994) 

[33] P.J. Forrester. J. Math. Phys. 35, 2539 (1994) 

[34] T. Papenbrock. Phys. Rev. A 67, 041601 (2003) 

[35] E. Strahov and Y.V. Fyodorov. arXive: math-ph/0210010| (2002) 

[36] C.W.J. Beenakker. Rev. Mod. Phys. 69, 731 (1997) 

[37] T.H. Baker and P.J. Forrester. J. Stat. Phys. 58, 1371 (1997) 

[38] J.P. Keating and N.C. Snaith. Commun. Math. Phys. 214, 57 (2000) 

24 



[39] CP. Hughes, J.P. Keating and N. O'Connell. Proc. R. Soc. A 456 2611 (2000) 
[40] J.P. Keating and N.C. Snaith. Commun. Math. Phys. 214, 91 (2000) 

[41] J.B. Conrey, D.W. Farmer, J.P. Keating, M.O. Rubinstein and N.C. Snaith. Commun. Math. Phys. 
237, 365 (2003) 

[42] CP. Hughes and Z. Rudnick. J. Phys. A 36, 2933 (2003) 

[43] J.P. Keating and N.C. Snaith. J. Phys. A 36, 2859 (2003) 

[44] P.J. Forrester. J. Phys. A 32, L159 (1999) 

[45] B. Rider. arXive |math.PR/0312043| (2003) 



25 



